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Abstract. Using Grassmann variables and an analogy with two dimensional elec- 
trostatics, we obtain the average eigenvalue distribution p{uj) of ensembles of N x N 
asymmetrically diluted Hopfield matrices in the limit — oo. We found that in the 
limit of strong dilution the distribution is uniform in a circle in the complex plane. 
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Random matrix theory has become an active field of research in mathematics and 
physics in the last decades. In physics, since the now classical work of Metha[|l| on the 
statistical description of the energy levels of atomic nuclei, random matrices have emerged 
as an important tool in the study of the localization transition[^ |^, quantum chaos|^, 
spin glasses!^, neural networks 0], and disordered systems in general. Most of the work 
deal with ensembles of hermitian or symmetric matrices whose individual properties are 
well known and can be exploited in more complex situations. In the last years the interest 
in asymmetric matrices has grown, motivated, besides its intrinsic mathematical value, 
by problems of dissipative quantum dynamics [§ and the dynamics of neural networks 0. 
In a recent paper, Sommers et.al.f^] calculated the average density of eigenvalues p{uj) 
of N X N random asymmetric matrices in the limit N —>■ oo, with elements Jij, given by 
a Gaussian distribution with zero mean and correlations 

N{{4))j = 1, iV((J,, J,,))j = T (1) 

They found that the eigenvalues are uniformly distributed inside an ellipse in the complex 
plane, whose semi axes depend on the degree of asymmetry of the ensemble r. Gener- 
alizing this result, Lehmann et. al. [TT| calculated the joint probability distribution of 
eigenvalues in Gaussian ensembles of real asymmetric matrices, recovering the elliptic law 
in the large N limit. 

In this paper we calculate the average eigenvalue distribution p(A) of an ensemble of 
asymmetrically diluted Hopfield matrices, whose elements are given by: 

j^^ = '-§i:m ^,j = i...iv, (2) 

where {^f « = l...N,fi = l...p} represents a set of p random patterns. The .^f are 
random independent variables that can take the values ±1 with the same probability and 
the Cij are random variables chosen accordingly to the following distribution 

P(q,)=75(%-1) + (1-7)5(q,). (3) 
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< 7 ^ 1 measures the degree of dilution of the matrices. 7 = 1 corresponds to symmetric 
Hopfield matrices whose eigenvalue distribution is known 

In order to obtain the distribution p{uj) we use an analogy with a two dimensional 
electrostatic problem introduced in ref.[|l^]. Let us define the Green function associated 
to the matrix J 

where uj = x + iy is a, complex variable, / the identity matrix and ((. . denotes an 
average over the random variables .^f . If Aj, z = 1, . . . , are the eigenvalues of J, then 

1 ^ 1 

Tr^^, = E (5) 
loj - J ^ uj - Xi ^ ' 

For large TV the sum can be approximated by an integral and the Green function becomes 

G{u) = I d'X-^ (6) 
J CO — X 

where p(A) is the density of eigenvalues in the plane. The last equation suggests an 
analogy with a two-dimensional classical electrostatics problem in which p(A) represents 



the density of charge in the plane. It can be demonstrated |[T0[] that an electrostatic 
potential $ exists, satisfying 

2ReG = -T^, -2ImG = -— (7) 

ox oy 

and which obeys Poisson's equation: 

= -Arrp (8) 

Thus, in order to determine p{uj) we may calculate the potential $. Using that det{AB) = 
detAdetB and detA^ = detA one can prove that a good definition for $ can be 

= -1/Ar((lndet ((/cj* - J^){Iuj - J)) ))j (9) 
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with uo* the complex conjugate of uo and the transpose of J. In what follows we will 
consider the case N oo and assume that in this limit the average and the In operations 
commute By using a Grassmannian representation [|13[ of the determinant and adding 



a matrix e 6ij, with e positive and infinitesimal in order to avoid zero eigenvalues, we get: 

roc / N \ 

exp [-iv$M] = ((y_^ n^^r^^O 



exp < - ^ Vi{^*^ik - Jlk){^^kj - Jkj)Vj \ ))j (10) 

After performing the average over the Qj and over the random patterns {^f }, we arrive 
at the following expression: 

exp [-N^{uj)] = 

P (jldmdvl^ exp{(e+ I u \^)Nq - aA^lnt + 07(0; + tu*)A^g - a7(l - 7)| u pA^g^} 
^/^ (X{d^.dXi^ exp |- X*Xi[^ + 07(1 - 7)5] + 

XiVi[<^7 - 07(1 - -f)ujq] + J2 V*iXi[ai - 07(1 - 7)0^*5] \ 

i i J 

X exp aATln 1 1 - -f\t J2 XlXi/N - 7t(E XlVi + E V:x^d/N+ 

K i i i 

7't(E ^Xl E mXj + E ^IXi E X*m)/N^ \ (11) 



where 



? = ^HvlVi (12) 

i 

1 . s 

t = ^7 ^ 13 



respectively. Next we define the following parameters: 

i=l i=l 



N N 

i=l 1=1 
^ N ^ N 

i=l i=l 

and introduce them into Eq. (|TT]) by using delta functions. After integrating over the 
Grassmann variables we get: 



exp [-N^{uj)] = 

— / dq dQ dz dZ dr dR dr* dR* 
X exp {gQ + zZ + rR + r*R* + [e+ \ oo |^ +a-f{uj + uj*)]q 
+a ln[l - 7(cu + uj*)q] + ln{RR* - ZQ) -z + a-f{r + r*) 
—0:7(1 — 7) [I uj pg^ + zq + ur*q + LJ*rq] 



+a In 



1 - 



'-y'^qz 



7(r + r*) 



+ 



9 ^ 



(15) 



1 — 'y{uj + uj*)q 1 — 7(co' + uj*)q ' 1 — 7(c<j + uj*)q 
In the large A^ limit this multiple integral can be evaluated by the saddle point method. Up 
to now the calculation is exact for arbitrary 7. Since the resulting saddle point equations 
are difficult to solve analytically, in this work we present the results for the strong dilution 
limit (7 <^ 1). Expanding the exponent in powers of 7 and keeping terms up to 0(7) we 
obtain, after some calculations: 



exp [-N<^{ijo)] oc 



+00 



(igexp — A^ In | g | — (e + | ui p)g + 0:7 | |^ — ln(l + Q;7g) 



(16) 



After a change of variables o" = 1/g we arrive at the following saddle-point equation: 



1 



X 



(17) 



cr^ cr + 07 (cr + a^Y (o" + 07)^ 
Expanding e in powers of a, the solution of the saddle-point equation in the limit e ^ 
is 0" = inside the circle + = 0:7. In this region G{uj) = uj*/ (0:7) (non analytic) and 
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= — 4/(«7). This implies that the density of eigenvalues is uniform inside a circle 
of radius ^^07 in the complex plane. Outside the circle the solution to Eq. ([T7|) becomes 
o = ^ — 0:7, the Green function is GiiS) = 1/uj (analytic), and the density p = 0. 
The density of eigenvalues in the whole complex plane is: 



p{uj) 



l/7ra7 if + y"^ < ^^^^ 
otherwise 



It is important to note that {{JijJji))j oc 7^ and consequently, in this limit of strong 
dilution, the matrix elements become effectively uncorrelated and we obtain a "circular 



law". This result has to be compared with the similar result of ref. [|10[ for the case r = 0, 
i.e., a Gaussian ensemble of completely asymmetric random matrices. It is expected that 
this circle deforms into an ellipse as the asymmetry parameter 7 increases and permits 
the appearance of random correlations between the patterns. 

In figures (1-3) we show the results of numerical diagonalization of sets of N x N 
matrices, with linear sizes ranging from 512 to 2048 for a = 0.25 and 7 = 0.01. The 
figures show the projection of the distribution of complex eigenvalues in the real axis. 
The full curves represent the analytical solution: 

p{x,y)dy 

2 

-(a7-2;2)V2^ \x\<^dn. (19) 



Tra'f 

We found that the numerical results present a peak at the origin that becomes smaller 
as the size N increases. Assuming that it is a finite size effect, and that the weight of 
the peak might be uniformly distributed on the whole support of the distribution, we 
renormalized the distributions. After renormalizing the numerical data the agreement 
with the analytic curves becomes very good as the size increases. Figure (4) shows the 
dependence of the peak at the origin with the system size. We have fitted the data at 



5 



the origin with an exponential function in 1/iV, Px{0) = aexp (b/N). The extrapolation 
to A?" — > oo coincides with an error of 10% with the analytic result at the origin. 

Concluding, in this paper we have obtained analytically the distribution of eingenval- 
ues of an ensemble of asymmetrically diluted Hopfield matrices in the limit of strong dilu- 
tion. The eingenvalucs are uniformly distributed inside a circle in the complex plane. Our 
results are supported by numerical diagonalization of the ensemble considered. Although 
we presented only the results for the strong dilution limit, the saddle point equations are 
valid for any amount of dilution. Work for solving the general case is in progress. 
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Caption for figures 

Figure 1 Projection of the eingenvalue distribution on the real axis. The full line corre- 
spond to the analytical results and the histogram to the numerical diagonalization 
performed with 7 — 0.01, a — 0.25 and N — 512 and averaged over 20 realizations 
of the matrices. 

Figure 2 Same as figure (1) with N — 1024 and averaged over 10 reahzations of the 
matrices. 

Figure 3 Same as previous figures with N = 2048 averaged over 8 realizations of the 
matrices. 

Figure 4 Finite size scaling of the peak at the origin in the complex plane In p(0) versus 
1/N. 
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